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Resumen 


MERCADO, J.R., GUIDO, P., SÁNCHEZ-SESMA, J. ée 
ÍÑIGUEZ, M. Fórmulas para el coeficiente de arrastre y la 
ecuación Navier-Stokes fraccional. Tecnología y Ciencias del 
Agua. Vol. V, núm. 2, marzo-abril de 2014, pp. 149-160. 


Se quiere encontrar la relación entre la ecuación de 
Navier-Stokes fraccional y las fórmulas para el coeficiente 
de arrastre, como las de Kármán-Schoenherr, Prandtl- 
Kármán, y Nikuradse. Los cambios de escala producen una 
renormalización para las ecuaciones de la capa límite, que 
contiene la hipótesis esencial de la delgadez de dicha capa, 
y da lugar a una descripción multifractal. Se obtiene una 
generalización del resultado experimental de Blasius para el 
factor de fricción. Si se reajustan las relaciones del número 
de rasgos del multifractal, se infieren las fórmulas, objeto de 
este estudio, y se las representa como un bi-multifractal, lo 
que permite un camino analítico para el número de Reynolds 
crítico y señala a la de Kármán-Schoenherr como la fórmula 
apropiada para el límite a la derecha de la subcapa viscosa. 
Los reajustes se traducen en matizar las aproximaciones de 
la relación entre los números de Euler y Reynolds, o bien 
en los decaimientos relativos del coeficiente de arrastre. 
Se aplican los resultados a la descripción de la capa límite 
turbulenta y a las interacciones entre corrientes y fondos (en 
ríos, desiertos y huracanes). 


Palabras clave: ecuación Navier-Stokes, coeficiente de arras- 
tre, multifractal, capa límite. 


Introducción 


El propósito principal del presente artículo 
es encontrar la relación entre las fórmulas de 
Kármán-Schoenherr, de Prandlt-Kármán y de 
Nikuradse para el coeficiente de arrastre o el 
factor de fricción, y la ecuación Navier-Stokes 
fraccional. 

En las referencias Mercado et al. (2013 y 
2012) se presenta la ecuación Navier-Stokes 
fraccional. El planteamiento esencial se basa en 
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afirmar que los esfuerzos viscosos producen un 
flujo dispersivo de momentum, que se describe a 
través de una ley de Darcy fraccional, y que la 
divergencia del flujo dispersivo coincide con el 
cambio temporal de momentum de acuerdo con 
la ley de Newton. 

Como en el caso clásico, la simplificación 
importante de la ecuación Navier-Stokes frac- 
cional, que conduce a las ecuaciones de la capa 
límite, estriba en su espesor relativamente 
delgado, lo que a su vez implica que la velocidad 
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principal es en la dirección longitudinal, 
con un gran gradiente de velocidad vertical, 
comparado con el longitudinal, mismo que 
lleva a la velocidad a cumplir la condición de 
no deslizamiento en el fondo del canal; y por 
el contrario, con los gradientes de presión leves 
en la dirección transversal vertical, comparados 
con los fuertes en la longitudinal (Landau y 
Lifshitz, 1987). 

Los cambios de escala se obtienen a partir 
de reconocer la longitud y velocidad carac- 
terísticas horizontales y macroscópicas, y 
buscar sus complementos verticales, que pre- 
serven la forma para la conservación de la 
masa y el momentum. 

Se introduce el potencial de la velocidad, 
por lo que tanto la componente de la velocidad 
aguas abajo como también la componente 
transversal vertical se dan por medio del 
mismo. Entonces, la ecuación de momentum se 
expresa en términos del potencial de velocidad 
aguas abajo. 

El esfuerzo friccional se halla calculando la 
derivada fraccional de la función potencial y se 
obtiene que el esfuerzo presenta la forma de una 
potencia del inverso del número de Reynolds 
indexado, y depende del índice de ocupación 
espacial; además, decrece con la distancia 
horizontal a una potencia también dependiente 
del grado de ocupación espacial. La constante 
de proporcionalidad contiene una potencia de 
la viscosidad y del valor de la curvatura de la 
función subpotencial en el origen. Se observa, 
además, que cuando el índice de ocupación 
espacial tiende a la unidad, se recupera la 
expresión clásica del esfuerzo friccional para la 
capa límite laminar. 

En la aproximación Falkner-Skan se esta- 
blece el equilibrio por medio de un triángulo 
dinámico conformado por las fuerzas viscosa, 
inercial y la del gradiente longitudinal de la 
presión, a través de la velocidad libre o exterior 
dependiente de una potencia de la coordenada 
longitudinal. En contraste, en la aproximación 
de Blasius, el gradiente de presión está ausente 
debido a la anulación de la potencia de la 
velocidad exterior. 


Se obtiene la fuerza de fricción en las dos 
aproximaciones mencionadas, para ello se 
integra el esfuerzo de fricción; se le expresa 
en forma adimensional, introduciendo el 
coeficiente de arrastre o drag. Resulta que este 
coeficiente también se representa como una 
potencia del número de Reynolds indexado, 
potencia determinada en forma acoplada 
tanto por el índice de ocupación espacial 
como por la potencia de la velocidad exterior. 
La potencia está conformada por el producto 
de dos factores; el primero es una función 
creciente, de manera que entre mayor sea el 
índice de ocupación espacial mayor será el 
exponente del número de Reynolds indexado, 
lo que se consigue con el movimiento viscoso; 
por el contrario, entre mayor sea la turbulencia 
menor será el índice de ocupación espacial 
y menor será el exponente del número de 
Reynolds indexado. El segundo factor es 
creciente para el exponente negativo, lo que 
trasluce un gradiente longitudinal adverso y es 
también creciente para el exponente mayor que 
la unidad que refleja un gradiente de presión 
favorable. 

Además, se observa que la fuerza de fric- 
ción puede describirse como una derivada 
fraccional de un multifractal, donde el orden de 
la derivada depende del índice de ocupación 
espacial, en tanto el espectro de dimensiones 
del multifractal depende, en forma acoplada, 
tanto de la potencia de la velocidad exterior 
como del índice de ocupación espacial. 

Pueden hallarse también el coeficiente de 
Chézy como una magnitud alternativa y la 
velocidad adimensional como razón entre la 
velocidad media y la velocidad de corte. 

La ecuación de Navier-Stokes clásica 
adolece de un parámetro que refleje el carácter 
fractal del movimiento del fluido. En tanto, su 
versión fraccional lo expresa en el orden de 
la derivada que aporta la fuerza de fricción 
viscosa. Si se quiere ver a qué corresponde ese 
orden, se puede estudiar la fuerza de fricción 
sobre un fondo plano extenso, a través de 
las ecuaciones de la capa límite. El resultado 
engloba una expresión experimental de Blasius, 
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que se interpreta como un multifractal. Si se 
mantiene esa interpretación multifractal para 
las otras fórmulas, como la de Prandtl-Kármán, 
se obtiene que las dimensiones fractales no 
superan el valor 1, y es más pequeño en la 
medida que el movimiento sea más turbulento. 

Por tanto, la afirmación es la siguiente: la 
ecuación Navier-Stokes fraccional aplicada a 
la capa límite, con un escalamiento que refleja 
la delgadez de dicha capa, y que se interpreta 
como un multifractal, produce las fórmulas de 
Blasius, Kármán-Schoenherr, Prandtl-Kármán 
y Nikuradse, como formas que establecen 
vínculos de dependencia de distintos grados 
entre el número de Euler y el número de 
Reynolds. Cada una de ellas puede describirse 
con una dimensión similar a la de Blasius, 
dimensión que se hace cada vez más pequeña a 
medida que sea mayor el número de Reynolds. 

En la sección “Ecuación Navier-Stokes frac- 
cional”, se describe la ecuación Navier-Stokes 
fraccional, ya considerada en Mercado et al. 
(2012). Se toman en cuenta las simplificaciones 
que introduce la capa límite. Se hallan el 
esfuerzo, la fuerza de fricción y el coeficiente 
de arrastre. Se obtiene una generalización del 
resultado experimental de Blasius. Se analiza 
su estructura multifractal. Luego se consideran 
las fórmulas de Kármán-Schoenherr, Prandtl- 
Kármán y Nikuradse. Se encuentra una ex- 
presión para el número de Reynolds crítico. 
Por último, se aplican a las interacciones entre 
fluidos y superficie como entre ríos y sus lechos, 
entre los vientos y la arena, y entre vientos y 
océano, como en los huracanes. 


Métodos 
Ecuación Navier-Stokes fraccional 


El movimiento del fluido se describe desde 
el punto de vista Euleriano, considerando un 
volumen de fluido limitado por una superficie 
de frontera, con su momentum por unidad de 
volumen dado por pv. Por su importancia, se 
considera primero la interacción por fricción 


interna. El gradiente fraccional se expresa por 
vi pv, donde p es la densidad de masa; v, la 
velocidad; f, el índice de ocupación espacial, 
y M es la medida de mezcla de las distintas 
direcciones espaciales. La difusividad del 
momentum es la a-viscosidad cinemática v,, 
así que el flujo de Darcy del momentum es 
q, La razón de cambio del momentum por 
unidad de tiempo es la divergencia negativa, 
o convergencia, del flujo de Darcy; o bien, se 
escoge M, tal que el flujo sea proporcional al 
negativo del laplaciano fraccional (Mercado et 
al., 2013), lo que se muestra en (1), con a. =1 + f: 


qp=v.Vipa, 
a/2 
=-Va (-A) pv (1) 


En seguida se toma en consideración la 
contribución de las variaciones de presión al 
cambio del momentum del fluido a través de la 
fuerza que aporta el gradiente de presión, de 
tal manera que la suma del esfuerzo de fricción 
viscoso y la presión hidrostática conforma 
el tensor T = v Vi, pv — pl y da lugar a la ley 
de deformación. A continuación se incorpora 
una fuerza potencial externa, por unidad de 
volumen, del tipo — Vp(. Después se introduce 
la hipótesis de la incompresibilidad del fluido. 
Se hace explícita la derivada material que 
componga la variación local con la advectiva. 
Pero la exigencia de objetividad requiere la 
invarianza bajo los cambios de coordenadas, 
por lo que se debe modificar el aporte 
advectivo y surge el término de la vorticidad. 
Por último, la contribución a la fuerza inercial 
de la vorticidad se escribe al lado derecho de 
la ecuación como v x rotv, y puede imaginarse 
como originado en una fuerza externa 
que dinamiza la evolución del campo de 
velocidades a través de su vorticidad, entrando 
en contradicción con la fuerza viscosa; en 
tanto, el otro término se interpreta como una 
restricción que a lo largo de las líneas de 
corriente contiene la ecuación de D. Bernoulli. 
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El coeficiente v, puede compararse con la 
viscosidad turbulenta de Boussinesq: 


Ly - vaca) Fvevaron=v| q(v00)+2 09] (2) 


Como ya se comentó, las ecuaciones de la 
capa límite se obtienen de la ecuación Navier- 
Stokes fraccional por las simplificaciones que 
se inducen a partir de la premisa de un espesor 
relativamente delgado. Ahora la ecuación de 
la capa límite bidimensional se le considera 
en su versión estacionaria O permanente, junto 
con la conservación de la masa en su forma de 
divergencia nula, tal como se muestra en la 
ecuación (3): 


ud u+vo,u=w,oju-0, (p/p), 0,u+ov=0 (3) 


Donde el gradiente de presión se da por 


d 
0, (p/p) -4-_U, U =u4(x,y) 


A la velocidad principal, aguas abajo, se le 


y>o 


da la forma potencial, representándola a tra- 
vés del potencial de la velocidad y(u,v) como 
U =0,Y, VD = 04. En la aproximación Blasius 
se establece un equilibrio entre dos: entre la 
fuerza viscosa y la fuerza inercial, al asumir 
como nulo el aporte del gradiente de presión. 
En la aproximación Falkner-Skan está presente 
el gradiente longitudinal de la presión como la 
tercera contribución al equilibrio, permitiendo 
que la velocidad libre o exterior dependa de 
la coordenada longitudinal. Surge entonces la 
función subpotencial g(8) como solución de 
la ecuación diferencial de la ecuación Blasius 
fraccional o Falkner-Skan fraccional (Landau y 
Lifshitz, 1987; White, 2006). 

Se calcula el esfuerzo por Toy = u,0%0,y(u,0), 


siendo u, = Ppv,. Con el coeficiente de 


Blasius  fraccional  B . Yet , donde 
B Ea ES 


A y 4) 


Cos y 


(3'(8)) , y el número de Reynolds 
E-0 


indexado Ri = ul*/y,, el esfuerzo puede 


escribirse como en (4): 


Ty =(8/0)8, 2 
(va) 


a 


Se calcula la fuerza de fricción por 
unidad de longitud binormal, o transversal- 
horizontal, como Fs =2 f, 7, da y con el 
esfuerzo representado como en (4), puede 


BY 1 
[o 20 
escribirse F¿=2 Loi B UN pu “lx “dx 
f a l/a | BJo 
(v,) 
(Mercado, 2010). En la aproximación 
Falkner-Skan se considera U =  Kx", 


con m % 0; se denota el coeficiente por 


E al Bx1U/m 1/a , 2 
ae = or iós ] Cog después 
de integrar se obtiene: 

2-P(1+1/m8) E 
Es =3PC 901) 65) 


Se adimensionaliza la fuerza de fricción por 
unidad de área de la superficie tangencial, por 
lo que surge el coeficiente de arrastre (6): 


B 
1 e (1+1/m8) 
e 


Al exponente del coeficiente de arrastre se 
le llamará el exponente Blasius (7): 


0(8.m)= E (1+1/m8) (7) 


Además, en la aproximación Blasius se 
considera m = 0; se integra y se obtiene el 
exponente Blasius: 


B 
0(8B,m=0)=-=— 8 
(Bam=0)- E (8) 
Se pueden reunir las dos en una sola 
expresión (9): 


Lp (1+1/m8), m0, 
0(p,m)= 


18 m=0, (9) 
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El exponente Blasius, en tanto dimensión, 
debe ser 0 >0, por lo que se requiere o bien que 
m > 0, lo que significa un gradiente de presión 
favorable, o bien m < -1/f, que se traduce en 
un gradiente de presión adverso. 

En particular, el estado dado por B => 0 
es equivalente al estado m = -1/f, y podría 
llamársele turbulento. 

Se observa también que es posible la 
representación O = Blo + 1/(1 + f)), porque para 
o > 0 existe un m > 0, dado por a 


op(1+8) 


que permite tal representación, misma que 
puede considerarse una generalización del 
resultado experimental de Blasius. Para m = 0 
existe la representación como fB/(1 + f), el cual 
contiene, entre otros, el resultado experimental 
de Blasius para $ = 1/3, y también el resultado 
para el régimen laminar 0 = 1/2, para f = 1 
(Mercado et al., 2013). Se aprecia que el valor 
de O > 0 puede llegar a ser bastante grande, 
a condición de tomar a m > 0 suficientemente 
cercano a cero. 

En la aproximación Falkner-Skan, para 
poder definir la variable de similaridad, 
se presenta la disyuntiva: om = 1,0 f = 1, 
porque son las dos opciones bajo las cuales 
los coeficientes de la ecuación, para la función 
subpotencial, se hacen independientes de x. 

En la primera opción, m = 1, se tiene 0(f,1) 
= 1, por lo que resulta que la variación de la 
presión con la velocidad es lineal como en el 
modelo de Hagen-Poiseulille para flujo laminar. 

Para la segunda opción, el exponente de 


Blasius es 0(1,m) = > dla 1 
m 


y el valor m es 
arbitrario, aunque, por supuesto, está excluido 
el valor m= 0; y además se requiere m > 0 o bien 
m<-—1/P, que en este caso se reduciría am<-1. 
La primera forma parte de los flujos inviscosos 
contra una elevación o cuña de ángulo 


a > 0; la segunda, de los flujos inviscosos, 
m+ 


a través de una depresión de ángulo T=s0, 


m+1 
en sentido horario (White, 2006). Además, 


se destaca que este exponente contiene el 
resultado experimental de Blasius con m = 2. 
Ahora la ley de variación de la presión con la 


velocidad es turbulenta, cercana a la fórmula 
de Chézy, la que se asume exactamente cuando 
se toma el valor m =-1. 

Es posible representar el factor de fricción 
por la relación que se muestra en (10) y descri- 
birlo como una transformación, por medio de 
derivadas fraccionales de orden y = f?/(1 + f) 


de un multifractal de dimensiones of = qe 

en otro análogo de dimensiones 6: mí Z P) 
fo=8C(p)D*(1 Ra 

9=8(1/mp(1+B)+1/(1+B)) (10) 


Blasius y un multifractal 


El término multifractal se usaba en general 
antes de 1980 y fue acuñado en turbulencia por 
Frisch (1995). Las propiedades multifractales 
son destacables cuando una distribución se 
muestra en conjuntos que exhiben la carac- 
terística de ser irregulares y esparcidos. Se trata 
de una medida que da la probabilidad de que 
un punto caiga en un conjunto determinado, 
pero con una distribución singular, de tal 
manera que no existe una densidad para 
describirla, y es entonces cuando el formalismo 
multifractal toma relevancia. Las propiedades 
multifractales constituyen el contexto apropia- 
do para describir las cualidades de una ley de 
escalamiento determinada. 

Se conforma una malla como una colección 
de cubos (C,), de lado 1/R, con 0 <1/R,< 1, 
denominando a 1/R, como la resolución. En el 
nivel microscópico, se considera una medida 
u, con la capacidad de estimar la probabilidad 
de que un punto determinado caiga en el cubo 
C, por u(C,) = (1/R ) bajo la condición natural 
de que la malla de cubos intersecte el soporte 
de la medida. Ahora se consideran todos los 
cubos en donde la medida sea del orden f, 
así que (1/R )' < u(C,) < (1/R*, siendo 1/R, 
pequeño, y e arbitrario, positivo y pequeño. 
Se le llaman rasgos al número de cubos en 
donde la medida es del orden de f (Falconer, 
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1990; Riedi y Scheuring, 1997; Meneveau y 
Sreenivasan, 1991): 


Na, (8)-(k:1(C,)=(1/R,)'] (11) 


La ley de escalamiento queda definida por 
la particularidad de que los rasgos N, (B) son 
del orden de 0(f) y, por tanto, obedecen a una 
ley de potencia, cuando la resolución se acerca 
a cero: 


Na (8) (Ry (12) 


Afirmación 1. El exponente de Blasius 
presenta una estructura multifractal. 

Las propiedades de un multifractal se pue- 
den analizar en el estimado que surge de la 
De Se 


debe recordar que el orden de la singularidad 


aproximación de Blasisus 0(8,m =0) = 


o dimensión local es f, que representa el índice 
de ocupación espacial o razón de la dimensión 
fractal a la dimensión topológica y el espectro 
multifractal de singularidades, como 0(f). La 
medida micro se escala como u = R;P, siendo la 
resolución 1/R, y R,elnúmero de Reynolds, y los 
rasgos como N, = RP), Se comprueba que 0(B) 
goza de la propiedad de la concavidad respecto 
de f3, lo que debido a su diferenciabilidad se 


puede inspeccionar a través de la segunda 
2 


derivada, porque (8) =0. La función de 
estructura es la transformada de Legendre 
del espectro que cambia por optimización 
la variable de la dimensión local f a la del 
orden de reiteración s u orden del momento. 
Siempre se tiene que t(s) es la transformada 
de Legendre del espectro aquí denotado por 


0(B) o sea sup(0(B)-sB) 
B 


cuando el espectro es diferenciable y se alcanza 


=1(s). En tanto que 


el supremo, la condición de optimalidad 


produce 59(6)=s. Se escribe 0(B) — sf = 1(s) 
P 


o —-sp="(s) y la condición de optimali- 


1+B 


dad arroja 7=5, así que invirtiendo se 


(1+p 
2 
tiene 1(s) = (Ys = 1) ; y entonces la condición de 
optimalidad pasa a ser (5) = -f, que tiene la 
s 
l OS J 

propiedad de la convexidad 1 = ——=3 

ds 2 E ] 

con la condición de normalización 1(1) = 0. 
¿Qué 
función de partición? 


>0 


magnitud física representa la 
Sis) = (1/RIO, 


2 


. Por 


NES: 2 2 B 
(gy y_ (2,0) , luego 5905 
B 

tanto, la función de partición representa al 


cuadrado del inverso del factor de fricción. 
Así que tomando como guía los fractales, 
se representa al factor de fricción por 


m5) 
fs 
In(R,) ! 


a1l/ R, como la resolución, bajo la condición 


donde se ve 


6 
4 =8Bg A o bien 8 = 


8B 
0<1 /R, <1, y ala cantidad —É como al número 


B 
de rasgos, por lo que se obtiene un multifractal 


de dimensión 0, que varía con f. Análogamente, 
e 2B; 
el coeficiente de arrastre —Y corresponde 


al número de rasgos de un multifractal de 
dimensión 0(f). También el coeficiente de 


A7n ar a 5) 


corresponde al número de 


Chézy satisface 


dond ST 
onde E sI [Bs 


rasgos de un multifractal de dimensión 8/2, y 
con resolución 1/R; 


(13) 


-8/2 
yu _ pl 
1 Js (2) 
Kármán-Schoenherr 


Los resultados experimentales del coeficiente 
de arrastre sobre una superficie plana presentan 
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una primera fase con exponente 1/2, para bajos 
números de Reynolds, resolución gruesa o alta; 
luego una región de transición con crecimiento 
ligero del coeficiente de arrastre, una meseta; 
luego un descenso que se aproxima por un 
exponente tipo 1/5, para resolución media, 
y después la representación de Kármán- 
Schoenherr (Rouse, 1946). 

La ley de Blasius se generaliza e interpreta 
como un multifractal, tomando al inverso del 
número de Reynolds como la resolución, y al 
espectro de singularidades por el exponente 
de Blasius. Esta ley puede explicar las fases 
de exponentes 1/2 y 1/5, descritas en el 
párrafo anterior, e incluso la meseta en el lí- 
mite cuando el exponente tiende a cero 9 = 0, 
lo que puede interpretarse como turbulencia 
totalmente desarrollada, en donde se elimina 
la participación de la viscosidad a través de 
anular el exponente. 

Se puede reconstruir el resultado de 


Kármán-Schoenherr por medio de debilitar 
1 


A,[C 
los rasgos desde B/ C, hasta Ke vo con dos 
constantes A y K. En congruencia con el méto- 
do para el resultado de Blasius, se considera 


1 
In[K EN 
ahora la dimensión 1-0="="___- 
Ink, 
in(B/C;) 
siendo O = $ si además es A = 4.13loge 
n e 
y KB = 1, se reconstruye la fórmula Kármán- 
Schoenherr o KS en (14): 
1 


== 4.1310g10(R.C) (14) 


15 


Por tanto, la fórmula KS corresponde 
también a una aproximación que establece una 
dependencia débil entre el número de Euler 
el de Reynolds, dado por E ——33 

ves B(R.) 
con la particularidad de que entre más intensa 
sea la turbulencia mayor será el valor del 


exponente en el número de Reynolds 1 — 6. Este 
exponente depende tanto de f3 como de m, con 
la restricción m > 1 requerida para que O < 1. 

Se inspecciona ahora su decaimiento y se le 
compara con la misma para la ley de Blasius. 
Se aprecia que al disminuir la dimensión 0, 
el decaimiento de C, relativo a sí mismo, es 
cada vez menor, porque si 6, < 8,, entonces 
El 
c; 


E: 
de ; por tanto, si se disminuye la 
0) 8; 


f 


dimensión 0, disminuye también la velocidad 
de decaimiento. Por otra parte, en el caso 


de Kármán-Schoenherr, se encuentra que 


Cs E 1 1 ; así se está cambiando 

Cs  [1+1/(24,[€,))R. 

la dimensión Oz Eb Oxs[Cy) == A la 
f 


cual quedaría definida de manera implícita. 
Se observa que como A,[C, > 0, se tiene que 
9,, < L, y que es cada vez más pequeño entre 
menor sea C, lo que se presenta cuando el 
número de Reynolds es suficientemente alto, 
y por tanto se obtiene una menor inclinación 
o pendiente. Se podría introducir un entero 


1,5 Como el mayor entero mayor o igual a 


eventualmente, con el crecimiento del número 
de Reynolds se produce una caída más leve con 
Kármán-Schoenherr que con la ley de Blasius. 
Por ejemplo, en un intervalo de cambio del 
coeficiente de arrastre disminuyendo de 0.010 
a 0.001, el número n,¿ cambia de 3 a 9, y O,. 
disminuye también de 1/4a1/10. 


Fórmula Prandtl-Kármán 


Por otra parte, y en congruencia con el 
desarrollo anterior, para la fórmula Prandtl- 
Kármán o PK se asume un multifractal con 
resolución 1/R, con rasgos proporcionales a 


A E 
Pffye” VÍ 


y de dimensión 1 — 6. Entonces, 
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con A =2 loge, D =0.8 y PB = 1, se recupera la 
fórmula Prandtl-Kármán: 


y =210310[/FR.)-08 (15) 


En el caso de Prandtl-Kármán, se encuentra 
el decaimiento relativo del coeficiente de 


—L--8 PK a. siendo la dimensión 
R 
ÑA e 
Opx (e ) = ER. APP También se 


¿(1=1/(24 €) 


observa que como A,/C f>0, se tiene que 


arrastre 


0 


px < 1, siempre que =1, Así, 0,7 S€ 


1 
hace cada vez más pequeño entre menor sea 
C, lo que se presenta cuando el número de 
Reynolds es suficientemente alto; por tanto, 
también se obtiene una menor inclinación o 
pendiente. Se introduce un entero n,, como 


1 
el mayor entero mayor o igual a 2A,/C+, lo 


que da 0, = s0px (Cy). Por tanto, 


5(1+1) 
eventualmente, con el crecimiento del número 
de Reynolds, se produce una caída más leve 
con Prandtl-Kármán que con la ley de Blasius. 
Por ejemplo, en un intervalo de cambio del 
coeficiente de arrastre, disminuyendo de 0.010 
a 0.001 el número n 


PK 


disminuye también de 2/7 a 1/10. 


cambia de 6 a 8, y Dor 


Fórmula de Nikuradse 


Ahora, en vez de los rasgos de la ley de Blasius 
1 


B/ Cs se les debilita por NC pen la idea es 


que si Cc, disminuye, con el incremento del 


número de Reynolds la cantidad NC + baja, pero 
1 


i : AC 
lo hace con mayor intensidad la cantidad e”; 


In|NC ye 
así, 1-08 = y nuevamente bajo la 
mR, 
In[B/C 
hipótesis O = e) y NB = 1, se obtiene: 
Ink, 
Alm(R,)=— (16) 
y 


, 


Su decaimiento relativo se da por a =-Oy = 
F e 


con By (c ,)=24 C en este caso, se aprecia 
que 6, puede ser mayor que la unidad 


cuando [Cy >1/24; sin embargo, en caso 


contrario también será menor que la unidad 
como en las tres fórmulas anteriores, pero 
permanece la propiedad de decaimiento de 0,, 
con el crecimiento del número de Reynolds. 
Por ejemplo, en el rango para el coeficiente 
de arrastre entre 0.01 y 0.001, 0, cambia, 
disminuyendo desde 0.358 hasta 0.113. 

Por tanto, las distintas fórmulas pueden 
verse como consecuencias de la ecuación 
Navier-Stokes fraccional, su reducción a la 
capa límite, su interpretación multifractal, y el 
reajuste de sus rasgos o de sus decaimientos 
relativos, lo que permite disminuir sus 
dimensiones con el crecimiento del número de 
Reynolds. Se puede entonces destacar al menos 
un bimultifractal tal, que para relativamente 
altos números de Reynolds los rasgos son pro- 

1 
porcionales a Ke We la resolución 1/R, y 
las dimensiones 1 — 6; en tanto que para 
bajos números de Reynolds, los rasgos son 
inversamente proporcionales al factor de 


oo 8B, : od 
fricción: —E, la resolución también 1/ R, y las 


B 
dimensiones 6: 


Nz = (17) 
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En especial, el cambio en la representación 
permite obtener un camino analítico para en- 
contrar un número de Reynolds crítico que 
señale el cambio de régimen de laminar a tur- 
bulento, o más específicamente, de la subcapa 
viscosa a la subcapa de mezcla o transición. 

Si el cambio de régimen se define a su vez por 
el cambio en la ley de la variación de la presión 
con la velocidad, se traduce en el cambio de la 
dimensión desde 9 = 1 a un valor inferior 6 < 
1 (Sommerfeld, 1950). Así, para números de 
Reynolds considerados relativamente bajos 
se tiene el comportamiento dado por la ley 
de Blasius, en tanto que para relativamente 
altos podría tomarse la fórmula de Kármán- 
Schoenherr, y el número de Reynolds que 
realiza la transición puede ser considerado 


como el Reynolds crítico. Así, en la fórmula KS 


se toma 0 = 1, se infiere que o =B;y conf 
= 4C, se obtiene f= 4/(A In B); en tanto de la 
fórmula de Blasius f=4/(B/R,), con 8 =1, 

por lo que el número de Reynolds crítico 
que señale la transición se puede estimar por 
R,. = B(A ln BY, o bien puede expresarse el 
coeficiente B en términos de A y R, hacien- 


do uso de la función de Lambert W y se obtie- 


ne B= Bio E 


44? w?(,/R.. /24) 
B = 64 y A = 4.13 loge = 4.13/In 10, se obtiene 
R,,= 3 561.1, que se valora del mismo orden de 


. En particular, para 


magnitud que el valor experimental estimado 
en 2000, y se recuerda que tanto el valor 
de B como el de A se determinan de modo 
experimental (Rouse, 1946). 

Por otra parte, en el caso de la fórmula 


de Nikuradse se combina 2/ AÍF =In4B/f 
con f = 4B/R, y se obtiene ¡Ro / AB =InR,, 


1 


1 1 
In = 
que se traduce en ,/R, 


JR.)  2AWB, 
cuya solución se expresa por la función de 
Lambert W, como Ro, = 4A4?BW?/(-1/2AWV/B; 
su evaluación numérica con B = 64 y A =4.13 / 


In10 produce R,, = 1.0749, el cual está en otro 
orden de magnitud, comparado con el valor 


2000. Análogamente, para la fórmula KP se 
obtiene, también a través de la función de 
Lambert, una estimación que produce 0.86709. 
En conclusión, se señala a la fórmula KS como la 
apropiada para el límite inferior de la subcapa 
de mezcla o transición. 


Resultados y discusión 


Se consideran tres aplicaciones. En la primera 
se quiere aportar en la descripción de la 
estructura de la capa límite turbulenta. En 
la segunda se estudia la interacción entre la 
rugosidad de la arena y la velocidad del viento. 
En la tercera se trata de ver la interacción entre 
velocidad del viento y rugosidad oceánica. 


Capa turbulenta 


Dentro de las aplicaciones se quiere aportar en 
la descripción de la estructura de la capa límite 
turbulenta. Se conocen cuatro zonas de la capa 
límite: una, la zona de pared, compuesta por 
dos subcapas: la viscosa y la de transición o 
mezcla. Luego la capa inercial con la turbulencia 
totalmente desarrollada y la estela. 

Para superficies lisas o regulares, la sub- 
capa viscosa puede extenderse hasta unas 
cinco unidades de pared o cinco longitudes 
características. La segunda, la de mezcla, hasta 
unas 30 longitudes características. La tercera, 
desde unos 30 hasta 400 unidades de pared. 
El resto corresponde a la estela. En tanto, para 
superficies rugosas, la longitud característica 
decrece, porque debe haber un desplazamiento 
hacia abajo debido a que la rugosidad 
incrementa las pérdidas de momentum (Clifford 
et al., 1993; Landau y Lifshitz, 1987; Levi, 1989). 

El multifractal para números de Reynolds 
relativamente bajos o intermedios se describe 


8B 
por ESE] y además la velocidad adi- 
Pp 


: u 8 
mensional es —=|=, por lo que se 
u. $ 
propone para la velocidad adimensional la 


rr) 
representación U,/ ¡Bo 2. Ahora, la 
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altura y se adimensionaliza con la longitud 


Pa A y v 5 
Característica O intrínseca 7% y se obtiene 


u, : 
—Y como el número de Reynolds asociado 
hd 

con la velocidad de corte, que se toma para 


la resolución del multifractal de dimensión 


8/2 
Ma E , que podría 


U. /,[B, Va 


B/2; entonces, 


escribirse U, = 
se recupera la aproximación de la ley de pared 


(y). Se observa que si O > 2, 


2 


ad es creciente po 8 > 0. Para la curva- 


. y 9 sl 
U, = y - La pendiente es y 7 4 


+ 


tura La, = ab - al 7 , entonces muestra 
y 


una primera subfase convexa si O = 2; e 
incluso para O = 4 se tiene el comportamiento 
parabólico. 

Una segunda subfase cóncava para O < 2 
puede ser representada como una ley tipo 
Blasius U = y!”, la cual es cóncava. Así que 
es posible describir la capa de mezcla desde 
el punto de inflexión hasta la capa totalmente 
turbulenta, y después también la estela o 
subcapa intermitente (Barenblatt et al., 1997): 


8B 
e In| —E 
u_ 1 (Uy? 0_ IN 3 8) 
u, gl vo)" 2 2MR,  2InR, 


Para la parte totalmente turbulenta se toma 
la fórmula de Nikuradse, con lo que se tiene 


a EPA así U, =V2AIn(y,), 
UN j 
la cuales lineal con el logaritmo, — u,= ADA. 

UN Y, 


5 cóncava. 


2 
creciente y sa JU J2A 
E (y) 


Parasuperficierugosa, setomanlarugosidad 


y la velocidad de corte como magnitudes para 
adimensionalizar la pendiente hidráulica J, 
y se garantiza la invarianza de forma de la 


U, es la velocidad de corte y g es la gravedad, 


que coincide en forma con la rugosidad de 
Charnock (1955) para la superficie oceánica 
(Mercado et al., 2012). 


Arena y viento 


En el transporte de arena por el viento 


fluctuante, la rugosidad también se describe 
2 
por y, =C; a con U,, la velocidad de corte; g 


la gravedad; y los valores del parámetro C,, se 
encuentran en el intervalo [0.02, 0.05] para flujo 
de arenas en túneles de viento y hasta 0.18 para 
experimentos de campo (Clifford et al., 1993). 
Aunque la fórmula es la misma, el parámetro 
C, tiene un rango de distintos valores que en la 
rugosidad de Charnock. 


Los huracanes 


U.y a 


Nuevamente = (u, ) o, esta rugosidad 


vo 

controla la variación vertical de la velocidad 
del viento (Powell et al., 2003), que lo equipara 
con el factor de fricción descrito por Nikuradse, 
suficientemente justificado para el régimen 
completamente rugoso, de tal manera que el 
rizado de la superficie marina producida por 
el viento clasifica a la superficie marina en el 
estado rugoso y permite aproximar el factor de 
fricción por la fórmula de Nikuradse. 

En el modelo de Barenblatt et al. (2005) se 
explica el fenómeno por la disminución de la 
energía turbulenta en la capa límite del fluido 
aire sobre la rugosidad oceánica, la cual se basa 
en la presencia de gotas de agua relativamente 
grandes O pesadas, lo que se traduce en la 
disminución del coeficiente de arrastre. 


Ríos y lecho 


La interacción entre el flujo y el lecho del río da 
lugar a estructuras coherentes en el lecho mis- 
mo, que presentan como una de sus caracterís- 
ticas cuantitativas la altura de la rugosidad, lo 


e ¡0N e 
que permite tomarla como frontera paa (u, Jo, 
vz 
ea u, 
O bien si ff = 1 tiene U, = a Y, 
va 
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El valor de la altura de la rugosidad puede 
hacerse corresponder con la condición de 
frontera para el perfil logarítmico de la capa 
inercial. De acuerdo con los resultados de 
Nikuradse, para superficie regular o lisa, el 
espesor de la subcapa viscosa es del orden de 
100 veces esa condición de frontera; en tanto, 
para superficie rugosa, el espesor de la subcapa 
es algo mayor (Rouse, 1946; Landau y Lifshitz, 
1987). Otra vez se adimensionaliza la altura 
y, de la rugosidad con la longitud intrínseca 


va 


, y se obtiene UY otra vez como número de 


» v 
Reynolds asociado con la velocidad de corte; 
pero debido a que la velocidad de corte para 
superficie rugosa es algo mayor que para la 
lisa, este número de Reynolds se incrementa. 


pb v v 
La longitud característica es 2=-2 


U. i/o 


su valor numérico con los datos citados en 
Clifford et al. (1993), A. Kirkbride, cap. 7., es 


6 
LAO aero? La longitud de la ru- 


./5/(1.000) : 


x* 


gosidad de Charnock (1955) es y, 7 ¡a 


5 


€ [0.015, 0.035]; pero ahora con valores de la 


constante en el intervalo [0.02, 0.05], e incluso 
con el valor 0.18. Sus valores se ubican en 


2 5 
(5/1 000) 0.02 1.0204 x 10 
E 0.05 2.551x10% | El 
0.18 9.1836x 10? 


espesor de la capa viscosa será del orden de 


1.0204x 107 1.0918x 107 
107x| 2.5551x107 |=| 2.7296x107? 
9.1836x 107 9.8265 x 107 


El valor reportado del espesor de la 
subcapa viscosa es de 2.1 x 107 m. Sin 
embargo, el límite inferior conocido para 
grava es de 2.0 x 10% m. Valor comparable 
con el segundo de los citados, que señala el 
coeficiente en C, = 0.05. La comparación entre 
la altura de la rugosidad y el espesor de la 
y, _ 1.0204x10” 
y. 107 
dentro del orden esperado. Sin embargo, para 


=102.4 está 


capa viscosa 


Y, 9.1836x10? 1 
S' 107 100 
puede desarrollarse debido a que la longitud 


, la subcapa viscosa no 


de la rugosidad excede, y con mucho, al espesor 
de esta subcapa. 


Conclusiones 


e  Alestudiar la interacción de un fluido con 
una superficie plana, a través de su capa 
límite, es posible establecer que la fuerza 
de fricción y, por tanto, el coeficiente de 
arrastre, se caracteriza por una potencia 
del inverso del número de Reynolds 
indexado, potencia que se ha denominado 
el exponente de Blasius, el cual presenta la 
dependencia de dos parámetros en forma 
de variables acopladas: uno manifiesta la 
viscosidad y el otro el gradiente de presión. 

+ De este exponente es posible establecer 
un estimado que responde a una estruc- 
tura multifractal, que se reduce al sufi- 
cientemente conocido para la subcapa 
viscosa (1/2), cuando el índice de ocupa- 
ción espacial se acerca a la unidad. 

+ Se incrementan los rasgos y se mantiene 
la estructura multifractal, con lo que es 
posible obtener las fórmulas de Kármán- 
Schoenherr, Prandtl-Kármán y Nikuradse 
para el factor de fricción para superficies 
planas y para tubos tanto lisos como 
rugosos. 

e Combinando la expresión multifractal con 
la velocidad adimensional, e introduciendo 
la altura de la rugosidad adimensional 
pueden describirse las cuatro subcapas de 
la capa límite turbulenta. 

e La 
un camino analítico para encontrar el 


descripción  bi-multifractal aporta 
número de Reynolds crítico, y solamente 
la fórmula de Kármán-Schoenherr pro- 
duce un estimado acorde con el orden de 
magnitud de resultados conocidos expe- 
rimentalmente. 
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